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INTRODUCTION 

Th»  molecular  machaniimi  of  both  !ir|«  itrain  plutic  deformation  and 
thoaa  encountered  in  internal  friction  in  flexible  chain  glauy  polymen  have 
remained  controvenial  at  beat.  While  many  ad-hoc  mechaniamt  of  ima|ined 
molecular  conformational  changes  (1-5)  have  been  advanced  and  have  proved  to 
be  very  uaeful  aa  kinetical  acalmj  laws  [l,4|,  there  hae  been  conaiderable  skep- 
ticiam  on  their  validity.  To  remove  thia  eontroveray  we  have  performed  detailed 
computer  aimulationa  of  auch  deformation!  on  bona-flde  molecular  atructural 
modela  of  |laaay  polypropylene  (PP)  obtained  earlier  by  Theodorou  and  Suter 
(T  L  S)  (6-aj,  and  on  biaphenol-A  polycarbonate  (PC)  which  we  have  newly 
determined.  Here  we  report  the  reaulta  of  theae  aimulationa. 


THE  MOLECULAR  STRUCTURE-MODELS 

The  molecular  atructure  model  ueed  for  {laeey  PP  waa  that  of  T  It  S  (4), 
for  flaeay  atactic  PP  of  a  degree  of  polymenaation  of  74,  at  a  denaity  of  0.I9Z 
j/cm1,  appropriate  for  a  temperature  of  233  K,  and  in  the  form  of  cubee  of 
14.2  A  on  the  tide,  aubject  to  periodic  boundary  conditiona. 

The  molecular  atructure  of  {laaay  PC  wae  obtained  by  a  similar  ipproath 
of  a  atatic  energy  minimiration  acheme  [9,10]  and  wai  for  a  PC  molecule  with 
a  degree  of  polymeritation  of  34,  containing  443  atomt,  in  the  form  of  a  cube 
of  18  44  A  on  the  aide,  having  periodic  boundary  condition!  and  at  a  denaity 
of  1.2  g/cm1,  appropriate  for  293  K. 

The  molecular  atructure  of  one  apecifle  configuration  of  PC  it  thown  in 
Pig.  1,  where  the  ateme  are  thown  reduced  to  about  1/3  of  their  tue  to  permit 
viewing  into  the  atructure. 


PLASTIC  STRAINING  QF  POLYPROPYLENE 

For  all  ataget  of  plutic  deformation  atomic  tite  atreaa  tenaora  were  calcu¬ 
lated  for  each  atom  by  the  generaliaatien  of  T  k  S  [4]  of  the  earlier  Born  at  i 
Huang  [11]  definition,  given  for  the  ith  atom  by: 


a;' 


where  o,  4  and  q  ue  carteaian  indices,  ■  and  ;  ue  indices  of  apecifle  atoms, 
r ij  is  the  radius  vector  connecting  stoma  ■  and  j,  Fff  and  T,]  we  the  4  and 
1  components  of  the  force  and  torque  vectors  acting  between  atoms  i  and 
while  f  *’  is  the  permutation  tensor,  and  Vt  is  the  van  der  Waals  volume  of 
atom  i.  Corresponding  to  the  atreases  atomic  site  strain  increments  were  also 
defined  based  on  the  increments  of  the  displacement  gradients  of  each  atom  in 
relation  to  its  immediate  neighborhood.  The  elemenu  of  this  strain  increment 
tensor  we  (for  infinitesimal  increments): 


At'* 


1  [3Au‘  4Au* 

2  dz*  +  dx* 


(*) 


Of  pwticulw  interest  were  the  trace  and  deviatoric  puts  of  the  atrsas  and  strain 
increment  tenaora,  giving  the  atomic  tite  negative  preuure  er,  von  Miaes  shew 
atreaa  r,  the  atomic  tile  dilatation  die,  and  the  work  equivalent  shew  strain 
increment  Aq,  defined  u: 

•  -  jrr(e)  :  r  as  jJ7l(e  -  el)*|,/*  P.«) 

Af=jrr(A<)  Aqr  =  |Tr|(i2i<  -  A<I)*]'/*J  (3,4) 

To  strain  the  initially  14.2A  simulation  cell  of  PP  under  a  constant  volume 
constraint,  tensile  strain  increments  of  2  x  10"J  and  two  associated  transverse 
compreuive  attain  increments  of  10~*  were  impoeed  on  the  cell.  All  atoms  in¬ 
side  the  cell  were  given  similar  affine  displacements  followed  by  re-minimisation 
of  the  energy  of  the  entire  structure. 


RESULTS  OP  PLASTIC  STRAINING  OF  PP  AND  DISCUSSION 

Figure  2a  shows  the  tensile  stress-strain  curve  of  a  single  configuration  de¬ 
formed  to  a  total  strain  of  0.2.  Because  of  the  very  small  site  of  the  cube,  it 
contains  a  substantial  level  of  atomic  ’stress  noise,*  ao  that  in  the  initial  state 
the  tensile  streu  is  not  tero.  Subsequent  deformation  steps  elicit  a  protracted 
antlutic  response  which  is  typical  of  glusy  structures  [12],  and  is  most  likely 
utociated  with  very  short  range  segment  relaxations.  Once  this  aneluticity 
is  completely  polarised  out,  the  structure  responds  in  a  series  of  clearly  delin¬ 
eated  cycles  of  elutic  loading  followed  by  plutic  ’collapse*  u  is  shown  at  b, 
d,  e,  f,  and  g.  Deformation  prior  to  a  plutic  collapse  exhibits  complete  elutic 
reversibility  while  reversal  of  strain  after  a  plutic  collapse  is  also  elutic  but  at 
the  newly  lowered  level  of  stress. 

When  the  responses  of  9  separate  configurations  ue  combined  as  in  Fig. 
2b,  the  level  of  stress  noise  is  very  substantially  reduced,  and  the  initial  tensile 
response  is  more  nearly  linearly  elutic.  Once  yielding  occurs  at  a  strain  level 
of  about  0.03  the  stress  strain  curve  continues  to  exhibit  individual  uncorre¬ 
lated  plutic  collapte  events.  Figure  2b  alto  shows  that  during  plutic  straining 
at  constant  volume  the  preuure  it  steadily  increuing  indicating  the  dilatant 
nature  of  the  plutic  deformation  in  the  glauy  PP. 

From  detailed  examinations  of  many  stereo  pairs  of  the  changing  molecular 
conformations  during  a  typical  plutic  collapse  event  what  could  be  gleaned 
wu  relatively  small  disorderly  displacements  of  ugments,  resulting  roughly  in 
a  quui-umform  small  plutic  strain  increment  over  the  entire  simulation  vol¬ 
ume.  Since  the  plutic  collapses  occur  at  a  constant  total  strain  the  elutic 
streu  drops,  Ao,  can  be  directly  related  to  plutic  strain  mcremenu,  Ac,,  by: 


Figure  3  shows  both  the  cumulative  distribution  and  the  frequency  distribu¬ 
tion  of  theu  plutic  strain  increments.  The  average  magnitude  of  these  strain 
increments,  when  converted  into  shew  strain  is  Aq‘  *  0.037(+0.035).  Theu 
are  nearly  an  order  of  magnitude  smaller  than  what  the  ewlier  ad-hoc  models 
had  considered,  where  it  had  usually  bun  auumed  that  relatively  luge  dis¬ 
placements  occur  among  relatively  small  groups  of  molecular  ugments.  The 
principal  finding  of  our  simulation  is  an  inverted  picture  in  which  a  plutic 
relaxation  event  consists  of  a  relatively  small  transformation  shew  strain  in  a 
relatively  luge  volume  element.  This  is  most  clearly  demonstrated  when  the 
simulation  results  art  compared  with  the  meuurements  of  the  activation  vol¬ 
umes  Av*  which  is  a  product  of  the  local  shew  strain  increment  Aq*  and  the 
site  of  the  local  volume  (I  in  which  a  thermally  activated  unit  plutic  relaxation 
event  occurs,  i.s., 

a»*  -  n  Aq*  *  *r  (fjp)  r-  (*) 

Table  I  gives  the  activation  volumu  meuured  by  Argon  and  Bessonov  [4] 
in  six  glauy  polymers.  If  the  transformation  ohtu  strains  of  our  simulation 
in  PP  art  taken  u  typical  for  most  glauy  polymers,  it  it  pouible  to  calculate 
estimates  of  the  situ  of  the  local  volumu  undergoing  a  unit  plutic  relaxation 
event.  Table  I  shows  that  theu  local  volumu  us  very  large  indeed,  being  typ¬ 
ically  40  -  40A  on  the  sids.  Clearly,  this  it  a  direct  consequence  of  the  relative 
inextentionality  of  backbone  bonds  and  inflexibility  of  bond  anglu,  requiring 
that  all  sheu  strains  develop  inside  this  lugs  volume  element  by  a  uriu  of 
succeuivs  bend  rotations  that  must  average  out  over  the  volume  to  a  value  no 
larger  than  the  overall  theu  strain  increment  in  the  volume. 

ROTATIONS  QF  PHENYLE.  2  RINGS  AND  CARBONATE  CRQUPS.1N.PC 

Streu  induced  rotations  of  phenylme  rings  and  perhaps  carbonate  groupe 
have  often  been  suspected  at  sourcu  of  the  4  internal  friction  peak  in  PC.  To 
explore  this  pouibility  most  phenylens  rings  in  the  simulation  structure  were 
incrementally  rotated,  followed  by  Axing  one  or  the  other  of  the  torsion  anglu 
on  either  side  of  the  rotated  phenylene  ring,  and  by  r*-minimising  the  system 
energy.  Such  a  ring  that  is  being  rotated  is  identified  in  Fig.  1  in  black,  In  both 
casu  the  total  energy  of  the  system  increased  stably  and  reversibly  to  a  peak 
value  where  upon  it  dropped  suddenly  to  a  lower  level.  Depending  on  whether 
the  torsion  angle  between  the  ring  and  the  cubonete  group  or  between  the  ring 
and  the  ieopropylidene  group  had  been  fixed,  two  generically  different  ruponsu 
were  found.  In  the  first  cau  the  drop  in  energy  wu  associated  with  a  complete 
■flip*  of  the  phenylene  ring  that  wu  being  rotated  while  in  the  second  cau  a 
deu  conformational  change  in  the  neighboring  eubenate  group  wu  found.  In 
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INTRODUCTION 

Th«  molacular  mechanisms  of  both  lar|<  attain  plaatie  deformation  and 
thoa*  oncountcrtd  in  internal  friction  in  flexible  chain  glassy  polymers  hava 
ramainad  controvarsial  at  bait.  Whila  many  ad-hoc  machaniima  of  imagined 
molacular  conformational  chengea  [1-5}  hava  been  advanced  and  hava  proved  to 
be  vary  uaaful  aa  kineticsl  scaling  lawa  'l,A],  thara  haa  baan  conaidarabla  akep- 
ticiem  on  their  validity.  To  ramova  thia  eonttovaray  wa  have  parformad  datailad 
computar  aimulationi  of  such  deformationi  on  bona-fide  molacular  atructural 
modala  of  glassy  polypropylana  (PP)  obtainad  earlier  by  Thaodorou  and  Sutar 
(T  k  S)  (6-8),  and  on  bisphanol-A  polycarbonata  (PC)  which  wa  hava  nawly 
determined.  Hera  wa  raport  tha  rasults  of  thasa  simulations. 


THE  MOLECULAR  STRUCTURE  MODELS 

Tha  molacular  structura  modal  usad  for  glassy  PP  wu  that  of  T  k  S  (6), 
for  glassy  atactic  PP  of  a  dagrea  of  polymantation  of  76,  at  a  dansity  of  0.692 
p/cm1,  sppropriata  for  a  tamparatura  of  233  K,  and  in  the  form  of  cubas  of 
18.2  A  on  tha  sida,  subjact  to  pariodic  boundary  conditions. 

Tha  molacular  structura  of  gluey  PC  wu  obtainad  by  a  similar  approach 
of  a  static  energy  minimization  schema  (9,10)  and  wu  for  a  PC  molecule  with 
a  degree  of  polymerization  of  31,  containing  183  atoms,  in  tha  form  of  a  cube 
of  18.11  A  on  tha  sida,  having  periodic  boundary  conditions  and  at  a  density 
of  1.2  p/cm1,  appropriate  far  29’/  K. 

The  molecular  structura  of  one  specific  configuration  of  PC  is  shown  in 
Fig.  1,  where  tha  atoms  are  shown  raducad  to  about  1/3  of  their  size  to  permit 
view'  ha  structura. 


HU.  AIXING  OF  POLYPROPYLENE 

Fc.  ,ei  of  "'•slic  deformation  atomic  site  stress  tensors  ware  calcu¬ 

lated  for  •  ate  a  generalization  of  T  k  S  [8j  of  tha  earlier  Born  and 
Huang  (111  Hefinit  for  tha  ifh  atom  by: 

•p-sigH*?"’?]  <■> 

where  a,0  and  q  are  cartesian  indicu,  i  and  /  art  indicu  of  spacific  atoms, 
rij  is  tha  radius  vector  connecting  atoms  1  and  j,  and  T(]  are  the  t  and 
7  components  of  tha  force  and  torque  vectors  ting  between  atoms  1  and  j, 
whila  e**1  is  tha  permutation  tensor,  and  V)  is  1  van  dar  Waals  volume  of 
atom  i.  Corresponding  to  tha  stresses  atomic  sits  strain  increments  wars  also 
defined  based  on  the  increments  of  tha  displacement  gradients  of  sach  atom  in 
relation  to  its  immediate  neighborhood.  The  elements  of  this  strain  increment 
tensor  are  (for  infinitesimal  increments): 


Ac*' 


1  dAu*  dAu* 

2  [1?"  +  ‘ 


(*) 


Of  puticular  inter*::  wars  the  traca  and  daviatoric  parts  of  the  stress  and  strain 
increment  tensors,  giving  tha  atomic  sits  negative  preuurs  t,  vor.  Mieos  shear 
streu  r,  tha  atomic  site  dilatation  At,  and  tha  work  aquivale.it  shear  strain 
increment  Aq,  defined  aa: 

9  *  jTr(s)  :  (3,1) 

Ar  =  j7V(At)  Aq  =  j77|(A<  -  A«I)‘|',,|  (5,6) 

To  strain  tha  initially  I8.2A  simulation  call  of  PP  under  a  constint  volums 
constraint,  tensile  strain  incramants  of  2  x  10*’  and  two  associated  transveru 
compressive  strain  increments  of  10"1  ware  imposed  on  tha  cell.  All  atoms  in¬ 
side  the  call  ware  given  similar  sfiine  displacements  followed  by  re-minimisation 
o(  the  energy  of  the  entire  structure. 


RESULTS  OF  PLASTIC  STRAINING  OF  PP  AND  DISCUSSION 

Figure  2a  shows  the  tansila  atrasi-itrain  curve  of  a  single  configuration  de¬ 
formed  to  a  total  strain  of  0.2.  Because  of  tha  vary  small  lisa  of  tha  cube,  it 
contains  a  substantial  level  of  atomic  ‘streu  noise,"  so  that  in  tha  initial  state 
tha  tensile  streu  is  not  zero.  Subsequent  deformation  steps  elicit  a  protracted 
analutic  response  which  is  typical  of  gluay  structures  (12),  and  is  moat  likely 
usociated  with  very  short  range  segment  relaxations.  Once  this  aneluticity 
is  completely  polarised  out,  the  structure  responds  in  a  aeries  of  clearly  delin¬ 
eated  cycles  of  elutic  loading  followed  by  plutic  “collapu"  u  is  shown  at  b, 
d,  e,  f,  and  g.  Deformation  prior  to  a  plutic  collapu  exhibits  complete  elutic 
reversibility  while  reversal  of  strain  after  a  plutic  collapu  is  also  elutic  but  at 
the  newly  lowered  level  of  stress. 

When  the  responses  of  9  aspirate  configurations  are  combined  u  in  Fig. 
2b,  the  level  of  stress  noise  is  very  substantially  reduced,  and  the  initial  tensile 
response  is  more  nearly  linearly  elutic.  Once  yielding  occurs  at  a  strain  level 
of  about  0.03  the  streu  strain  curve  continues  to  exhibit  individual  uncorre¬ 
lated  plutic  collapu  events.  Figure  2b  also  shows  that  during  plutic  straining 
at  constant  volume  the  preuure  is  steadily  increuing  indicating  the  dilatant 
nature  of  the  plutic  deformation  in  the  glauy  PP. 

From  detailed  examinations  of  many  stereo  pairs  of  the  changing  molecular 
conformations  during  a  typical  plutic  collapu  event  what  could  be  gleaned 
wu  relatively  small  disorderly  displacements  of  ugmenta,  ruulting  roughly  in 
a  quui-umform  small  plutic  strain  increment  over  the  entire  simulation  vol¬ 
ume.  Since  the  plutic  collapses  occur  at  a  constant  total  strain  the  elutic 
streu  drops,  A  o,  can  be  directly  related  to  plutic  strain  increments,  Ar,,  by: 


Figure  3  shows  both  the  cumulative  distribution  and  the  frequency  distribu¬ 
tion  of  theu  plutic  strain  increments.  The  average  magnitude  of  theu  strain 
increments,  when  converted  into  shear  strain  is  Aq*  *  0.037(+0.033).  Theu 
are  nearly  an  order  of  magnitude  smaller  than  what  the  earlier  ad-hoc  models 
had  considered,  where  it  hed  usually  bun  anumed  that  relatively  large  dis¬ 
placements  occur  among  relatively  small  groups  of  molaculu  ugmenta.  The 
principal-finding  of  our  simulation  is  an  inverted  picture  in  which  a  plutic 
relaxation  event  consists  of  a  relatively  small  transformation  shear  strain  in  a 
relatively  large  volume  element.  This  is  most  clearly  demonstrated  when  the 
simulation  results  are  compared  with  the  meuurementa  of  the  activation  vol¬ 
umes  A v*  which  is  a  product  of  the  local  sheu  strain  increment  Aq'  and  the 
site  of  the  local  volume  fl  in  which  a  thermally  activated  unit  plutic  relaxation 
event  occurs,  i.e., 


Av"-QA,'-*r(^i)>r.  (8) 

Table  I  gives  the  activation  volumu  meuured  by  Argon  and  Bessonov  (1] 
in  six  glauy  polymers.  If  the  transformation  shear  strains  of  our  simulation 
in  PP  are  taken  u  typical  for  most  glauy  polymers,  it  is  pouible  to  calculate 
estlmatu  of  the  sizes  of  the  local  volumu  undergoing  a  unit  plutic  relaxation 
event.  Table  I  shows  that  theu  local  volumu  are  very  large  indeed,  being  typ¬ 
ically  10  -  60A  on  the  side.  Clearly,  this  is  a  direct  consequence  of  the  relative 
inextensionality  of  backbone  bonds  and  inflexibility  of  bond  anglu,  requiring 
that  all  shear  strains  develop  inside  this  large  volume  element  by  a  uriu  of 
succeui-e  bond  rotations  that  mutt  average  out  over  the  volume  to  a  value  no 
larger  than  the  overall  shear  strain  increment  in  the  volums. 


ROTATIONS  QF  PHENYLENE  RINGS  AND  CARBONATE.gROUI?S.ltt£C 
Streu  induced  rotations  ol  phenylene  rings  and  perhaps  carbonate  groups 
have  often  bun  suspected  u  sourcu  of  the  f  internal  friction  peak  in  PC.  To 
explore  this  pouibility  most  phenylene  rings  in  the  simulation  structure  were 
incrementally  rotated,  followed  by  Axing  one  or  the  other  of  the  torsion  anglu 
on  either  side  of  the  rotated  phenylene  ring,  and  by  re-minimising  the  system 
energy.  Such  a  ring  that  is  being  rotated  is  identified  In  Fig.  I  in  black.  In  both 
casu  the  total  energy  of  the  system  inereeeed  stably  and  reversibly  to  a  peak 
value  where  upon  it  dropped  suddenly  to  a  lower  level.  Depending  on  whether 
the  torsion  angle  betwun  the  ring  and  the  carbonate  group  or  betwun  the  ring 
end  the  isopropylidene  group  had  been  fixed,  two  genetically  different  ruponeu 
were  found.  In  the  first  cau  the  drop  in  energy  wu  associated  with  a  complete 
•dip'  of  the  phenylene  ring  that  wu  being  rotated  while  in  the  second  cau  a 
dear  conformational  change  in  the  neighboring  carbonate  group  wu  found.  In 
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